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On Which Length Scales Can Temperature Exist in Quantum Systems? 
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We consider a regular chain of elementary quantum systems with nearest neighbor interac- 
tions and assume that the total system is in a canonical state with temperature T. We analyze 
under what condition the state factors into a product of canonical density matrices with respect 
to groups of n subsystems each, and when these groups have the same temperature T. While 
in classical mechanics the validity of this procedure only depends on the size of the groups n, 
in quantum mechanics the minimum group size n m i n also depends on the temperature T ! As 
examples, we apply our analysis to different types of Heisenberg spin chains. 
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Thermodynamics is among the most successfully 
and extensively applied theoretical concepts in physics. 
Notwithstanding, the various limits of its applicability 
are not fully understood. 2 ' 7,8 

Of particular interest is its microscopic limit. Down to 
which length scales can its standard concepts meaning- 
fully be defined and employed? 

Besides its general importance, this question has 
become increasingly relevant recently since amazing 
progress in the synthesis and processing of materials with 
structures on nanometer length scales has created a de- 
mand for better understanding of thermal properties of 
nanoscale devices, individual nanostructures and nanos- 
tructured materials. 1 Experimental techniques have im- 
proved to such an extent that the measurement of ther- 
modynamic quantities like temperature with a spatial 
resolution on the nanometer scale seems within reach. 3 ' 4 

To provide a basis for the interpretation of present day 
and future experiments in nanoscale physics and tech- 
nology and to obtain a better understanding of the lim- 
its of thermodynamics, it is thus indispensable to clarify 
the applicability of thermodynamical concepts on small 
length scales starting from the most fundamental the- 
ory at hand, i. e. quantum mechanics. In this context, 
one question appears to be particularly important and 
interesting: Can temperature be meaningfully defined on 
nanometer length scales? 

The existence of thermodynamical quantities, i. e. the 
existence of the thermodynamical limit strongly depends 
on the correlations between the considered parts of a 
system. 6 With increasing size, the volume of a region 
in space grows faster than its surface. Thus effective in- 
teractions between two regions, provided they are short 
ranged, become less relevant as the sizes of the regions 
increase. 5 ' 9 The scaling of interactions between parts of 
a system compared to the energy contained in the parts 
themselves thus sets a minimal length scale on which cor- 
relations are still small enough to permit the definition 
of local temperatures. It is the aim of this letter to study 
this connection quantitatively 



For this purpose, the applied definition of temperature 
becomes crucial. We define local temperature to exist if 
the respective part of the system is in a canonical state. 
Besides its motivation from statistical mechanics, there 
are also rather practical arguments for this definition. 

The canonical distribution is an exponentially decay- 
ing function of energy characterized by one single param- 
eter. 11 This implies that there is a one-to-one mapping 
between temperature and the expectation values of ob- 
servables, by which temperature is usually measured. 15 

Furthermore, the product of this distribution times 
the density of states, which is typically a strongly grow- 
ing function of energy 10 for large systems, forms a pro- 
nounced peak and thus physical quantities like energy 
have "sharp" values. 

Based on the above arguments and noting that a 
quantum description becomes imperative at nanoscopic 
scales, the following approach appears to be reasonable: 
Consider a large homogeneous quantum system, brought 
into a thermal state via interaction with its environment, 
divide this system into subgroups and analyze for what 
subgroup-size the concept of temperature is still applica- 
ble. 

Recently, spin chains have been subject of extensive 
studies in condensed matter physics. We therefore study 
various types of Heisenberg chains with respect to our 
present purpose. 

We consider a homogeneous (i.e. translation invariant) 
chain of elementary quantum subsystems with nearest 
neighbor interactions. The Hamiltonian of our system is 
thus of the form 

H = J2 H i + hi+u (1) 

i 

where the index i labels the elementary subsystems. Hi 
is the Hamiltonian of subsystem i and the interac- 

tion between subsystem i and i + 1. We assume periodic 
boundary conditions. 

We now form Nq groups of n subsystems each (index 
i -> (/x - l)n + j; fi = 1, . . . , N G ; j = 1, . . . , n) and split 
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this Hamiltonian into two parts, 

H = H + I, (2) 

where H is the sum of the Hamiltonians of the isolated 
groups, 



N G 



*■ fj,n,iJ,n-\-l ) 



with 



(3) 



Kfj. ~ ^ -g w ( M -l) + j + I n (n-l)+j, n (n-l)+j+l 

i=i 

and / contains the interaction terms of each group with 
its neighbor group, 



N a 

I ^ ^ Ifj,n,fj.n-\-l' 



(4) 



The eigenstates of the group Hamiltonian H are prod- 
ucts of group eigenstates, 

N G 

H Q \a) = E a \a) with \a) = JJ |a„), (5) 

n=i 

where (H^ - I^nn+i) |o^> = ^ is tne energy 

of one subgroup only and E a = Y^flZi 

We assume that the total system is in a thermal state, 



-0H 



(6) 



where Z is the partition sum and [3 — (fcsT) -1 the in- 
verse temperature with Boltzmann's constant ks and 
temperature T. 

For our purpose, we need to know the representation 
of p in the basis formed by the product states \a). The 
diagonal elements (a\p\a) are the expectation values of p 
in the states la), 



(a\p\a) = / w a (E) —— dE, 

J Eq 



(7) 



where w a (E) is the probability density of the occurrence 
of the energy eigenvalue E in the state \a). E is the 
energy of the ground state and E\ the upper limit of 
the spectrum. For systems with an energy spectrum that 
does not have an upper bound, the limit E\ — > oo should 
be taken. 

One can show that a quantum central limit theorem 
exists for the present model 12 ' 14 and that w a (E) thus 
converges to a Gaussian normal distribution in the limit 
of infinite number of groups Ng, 

lim w a (E) = * — exp (- {E ~ E " ; ~ £o)2 ) , (8) 

where the quantities e a and A a are defined by 

e a = (a\H\a)-E a (9) 

A 2 ee (a\H 2 \a) - (a\H\a) 2 . (10) 

e a is the difference between the energy expectation value 
of the distribution w a (E) and the energy E a , while A 2 , is 
the variance of the energy E for the distribution w a (E). 



Note that e a has a classical counterpart while A 2 is 
purely quantum mechanical. It appears because the com- 
mutator [H,Hq] is nonzero, and the distribution w a {E) 
therefore has nonzero width. 

The rigorous proof of eq. (8) is given in 12 and based on 
the following two assumptions: The energy of each group 
Hfj, as defined in eq. (3) is bounded, i. e. 

(x\n,\ x )<c (ii) 

for all normalized states |%) and some constant C, and 

(a\H 2 \a) - (a\H\a) 2 > N G C (12) 

for some constant C > 0. 

If conditions (11) and (12) are met, eq. (7) can be 
computed for Ng S> l: 13 

1 / 2 A 2 
(a\f>\a) = — expl~/3y a + !-^ 



rfc 



E Q -y a + /3A 2 a 
V2A a 



erfc 



E 1 -y a + (3 Aj 
V2A a 



2 \ 1 



(13) 



where y a — E a + e a and erfc(x) is the conjugate Gaus- 
sian error function. 16 The second error function in (13) 
only appears if the energy is bounded and the integra- 
tion extends from the energy of the ground state Eq to 
the upper limit of the spectrum E\ . 

Note that y a is a sum of Nq terms and that A a fulfills 
cq. (12). The arguments of the conjugate error functions 
thus grow proportional to ^Nq or stronger. If these ar- 
guments divided by V ' Nq are finite (different from zero), 
the asymptotic expansion of the error function 16 may 
thus be used for Nq 3> 1: erfc(x) — > exp (—a; 2 ) /(y/nx) 
for x — > oo and erfc(x) — > 2 + (exp (— x 2 )) /{yfnx) for 
x — » — oo. Inserting this approximation into eq. (13) and 
using E < y a < Ex shows that the second conjugate 
error function, which contains the upper limit of the en- 
ergy spectrum, can always be neglected compared to the 
first, which contains the ground state energy. 

The off diagonal elements (a\p\b) vanish for 
\E a — Eb\ > A a + Ab because the overlap of the two distri- 
butions of conditional probabilities becomes negligible. 
For \E a — Eb\ < A a + A&, the transformation involves 
an integral over frequencies and thus these terms are 
significantly smaller than the entries on the diagonal. 

We now test under what conditions the density matrix 
p may be approximated by a product of canonical den- 
sity matrices with temperature f3\ oc for each subgroup 
p = 1,2,..., Ng. Since the trace of a matrix is invari- 
ant under basis transformations, it is sufficient to verify 
the correct energy dependence of the product density 
matrix. If we assume periodic boundary conditions, all 
reduced density matrices are equal and their product is 
of the form (a\p\a) oc cxp(— 0\ oc E a ). We thus have to 
verify whether the logarithm of rhs of eq. (13) is a linear 
function of the energy E a , 



In {(a\p\a)) w -/3\ oc E a + c, 



(14) 



where j3\ oc and c are constants. 

We exclude negative temperatures ((3 > 0). Eq. (14) 
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can only be true for 

E a +e 



Eq 



>/3- 



(15) 



as can be seen from eqs. (13) and the asymptotic ex- 
pansion of the error function. To satisfy (14), e a and A 2 
furthermore have to be of the following form: 



2 a 



CiE a + c 2 



(16) 



where ci and c-i are constants. Note that e a and A 2 need 
not be functions of E a and therefore in general cannot 
be expanded in a Taylor series. 

To ensure that the density matrix of each subgroup \x 
is approximately canonical, one needs to satisfy (16) for 
each subgroup \i separately; 



' 1 + £ f +^ (A 2 , + A 2 ) + ^A 2 



(17) 



where 

v=u- 



and 



a" 



£ p = (o|/ M „, Mn+ i|a) with 

= (a|H 2 |a) - (a|W» 2 

_ x <a|W v _iW w +W w W 1 ,_i |a) - 2<o|W v _i|o) (a|W„|a). 
Temperature becomes intensive, if the constant ci van- 
ishes, 



Icil « 1 



Aoc - P. 



(18) 



If this was not the case, temperature would not be in- 
tensive, although it might exist locally. 

It is sufficient to satisfy conditions (15) and (17) for 
an adequate energy range i? m i n < < E max only. For 
systems composed of a large number of subsystems, the 
density of states is typically a rapidly growing function of 
energy. 10 If the total system is in a thermal state, occupa- 
tion probabilities decay exponentially with energy. The 
product of these two functions is thus sharply peaked at 
the expectation value of the energy E of the total sys- 
tem E + E =Tr(Hp), with E being the ground state 
energy. Therefore a pertinent and "safe" choice for E m i n 
and .E max is 



E n 

Em 



max 



mm 



1 E 

MJmin ' a No 
E 



Eg \ 
No j 



max ' ^ No N G ) 



(19) 



where a>l and E will in general depend on the global 
temperature. In eq. (19), [E^] min and [-E M ] max denote the 
minimal and maximal values E^ can take on. 

For a model obeying eqs. (11) and (12), the two con- 
ditions (15) and (17), which constitute the general result 
of this letter, must both be satisfied. These fundamental 
criteria will now be applied to a concrete example. 

We now consider the Heisenbcrg-model of a spin chain 
in a transverse field . The Hamiltonian reads 17 ' 18 

H l = Bol ■ h, i+1 = J {a* a*at +l ) 

(20) 

where of, of and of are the Pauli matrices. B is the 
magnetic field and J the coupling parameter. We will 
always assume B > 0. 

We partition the chain into Nq groups of n adja- 
cent spins, as considered above, and numerically an- 
alyze eqs. (15) and (17) by exact diagonalization of 



the groups. In doing so, we apply the approximation 
A 2 « A 2 . Since A 2 = E^tl-M^-iL + Uu-i\a) - 
2{a\I v -i\a){a\I v \a), only o z ® o z terms contribute and, 
among those, only terms which are products of one diag- 
onal element and one off-diagonal element of the two a z . 
A 2 thus contains much more terms, which, in addition, 
are all positive. 

The conditions for the central limit theorem are met 
for almost all states \a) of the present model: Condition 
(11) is fulfilled because the Hamiltonian of a single spin 
has finite dimension. The numerics shows that for almost 
all states, A 2 is a finite positive number, which implies 
that condition (12) is satisfied. 

For single spins, a local temperature can always be 
defined by assigning a Boltzmann factor to the occupa- 
tion probabilities of the upper and lower level. This local 
temperature, however, is not equal to the global one. 

For groups of more than one spin, neither nor A 2 
can be approximated by linear functions of E^. Therefore 
local temperature can only exists if the rhs of eq. (17) 
is a constant. As a direct consequence, temperature is 
intensive, (3\ oc — /3. 

Eq. (15) can be checked directly. To check eq. (17), we 
need to quantify what "approximately constant" means. 
We thus use 



[A 2 1 - [A 2 1 



- \E U ] ■ 

l pJmin 



(21) 



where [x] max and [x] m ;„ denote the maximal and minimal 
value x takes on in all states \a). Condition (21) implies, 
that (17) holds. 

By numerical evaluation of eqs. (15) and (21), one can 
calculate a minimal group size n m ; n > 2 for each temper- 
ature T, where we consider (21) to be satisfied if the lhs 
is 100 times smaller than the rhs. This corresponds to a 
tolerable deviation from the canonical distribution of 1%. 
Figure 1 shows function of T for antiferromag- 

netic coupling with strength J = B. Local temperature 
can exist in the shaded region. As mentioned above, local 
temperature can always be defined for single spins. For 
the present coupling strength, however, it is not intensive 
unless T — > oo. 

Figure 2 shows function of the temperature 

T for antifcrromagnctic couplings of different strength. 
Diamonds correspond to J = 0.1B, stars to J = IB and 
squares to J = 105. The regions, where local tempera- 
tures exist are as indicated in figure 1. For single spins, 
local temperature always exists, but is only intensive for 
J = 0.1B and T > 2. 55 

In both plots, only the one of eqs. (15) and (17), which 
sets the stronger bound enters. For the present parameter 
range, this is exclusively eq. (17). The results for the 
ferromagnetic case are approximately the same and are 
therfore not given separately here. 

In summary, we have considered a homogeneous chain 
of quantum systems with nearest neighbor interactions. 
We have partitioned the chain into identical groups of 
n adjoining subsystems each. Taking the number of such 
groups to be very large and assuming the total system to 
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Fig. 1. 



as a function of the temperature T for J = B > 0. 



T is given in units of B, ks = 1 and a = 10. a is defined in 
eq. (19). Local temperature exists in the shaded region. It also 
exists for single spins but is never intensive there. 



10 1 10 2 10 3 



T/B 



Fig. 2. n m ; n n m i n as a function of the temperature T for J = 
0.1B (diamonds), J = B (stars) and J = 10B (squares), J > 
and B > 0. T is given in units of B, kg = 1 and a = 10. a is 
defined in cq. (19). For each value of J, local temperature exists 
for n > n m ; n . For single spins, local temperature always exists, 
but is only intensive for J = 0.1B and T > 2.5B. 



be in a thermal state with temperature T we have found 
conditions (eqs. (15) and (17)), which ensure that each 
group is approximately in a thermal state. Furthermore, 
we have determined when the isolated groups have the 
same temperature T, i.e. when temperature is intensive. 

The result shows that, in the quantum regime, these 
conditions depend on the temperature T, contrary to the 
classical case. The characteristics of the temperature de- 
pendence are determined by the width A a of the distri- 
bution of the total energy eigenvalues in a product state 
and its dependence on the group energies E a . The low 
temperature behavior, in particular, is related to the fact 
that A a has a nonzero minimal value. This fact does not 
only appear in spin chains but is a general feature of 
quantum systems composed of interacting particles or 
subsystems. The commutator [H, H ] is nonzero and the 
ground state of the total system is energetically lower 
than the lowest product state, therefore A a is nonzero, 
even at zero temperature. 19 



For the models we consider here, the off diagonal ele- 
ments of the density operator in the product basis, {a\p\b) 
(a 7^ b), are significantly smaller than the diagonal ones, 
(a\p\a). Our general result, conditions (15) and (17), thus 
states that the density matrix p "approximately" factor- 
izes with respect to the considered partition. This implies 
that the state p is not entangled with respect to this 
partition, at least within the chosen accuracy. It would 
therefore be interesting to see how our result relates to 
the scaling of entanglement in many particle systems. 20 

Unfortunately, our approach only applies to nonzero 
temperatures. The underlying central limit theorem 12, 14 
is about the weak convergence of the distribution of en- 
ergy eigenvalues. Weak convergence means that only in- 
tegrals over energy intervals of nonzero length do con- 
verge. We thus cannot make statements about a system 
in its ground state let alone about the entanglement in 
that state. 

We have then applied the general method to several 
types of Heisenberg spin chains. For concrete models, the 
conditions (15) and (17) determine a minimal group size 
and thus a minimal length scale on which temperature 
may be defined according to the temperature concept we 
adopt. Grains of size below this length scale are no more 
in a thermal state. Thus temperature measurements with 
a higher resolution should no longer be interpreted in a 
standard way. 

The length scales, calculated in this paper, should also 
constrain the way one can meaningfully define tempera- 
ture profiles in non-equilibrium scenarios. 21 



9 
10 

11 

12 

13 

14 
15 
16 

17 
18 
19 

20 

21 



D. Cahill et al.: J. Appl. Phys. 93 (2003) 793. 

A.E. Allahverdyan and Th.M. Nicuwcnhuizcn: Phys. Rev. Lett. 

85, (2000) 1799. 

Y. Gao and Y. Bando: Nature 415 (2002) 599. 
H. Pothicr ct al.: Phys. Rev. Lett. 79 (1997) 3490. 

D. Ruelle: Statistical Mechanics W.A. Benjamin Inc., New 
York (1969). 

K. Saito, S. Takesue and S. Miyashita: J. Phys. Soc. Jpn. 65 
(1996) 1243. 

T.L. Hill: Nano Lett. 1 (2001) 273. 

A.K. Rajagopal, C.S. Pande and S. Abe: cond-mat/0403738. 
M. Hartmann et al.: Europhys. Lett. 65 (2004) 613. 
R.C. Tolman: The Principles of Statistical Mechanics Oxford 
Univ. Press, London (1967). 

R. Kubo, M. Toda and N. Hashitsume: Statistical Physics II 
Springer, Berlin (1985). 

M. Hartmann, G. Mahler and O. Hess: Lett. Math. Phys. 68 
(2004) 103, math-ph/0312045. 

M. Hartmann, G. Mahler and O. Hess: Phys. Rev. Lett. 93 
(2004) 080402, quant-ph/0312214. 

M. Hartmann, G. Mahler and O. Hess: cond-mat/0406100. 

M. Hartmann and G. Mahler: cond-mat/0410526. 

M. Abramowitz and I. Stcgun: Handbook of Mathematical 

Functions Dover, New York (1970) 9th ed. 

J. van Vleck: Rev. Mod. Phys. 17 (1945) 27. 

E. Manousakis: Rev. Mod. Phys. 63 (1991) 1. 

A.N. Jordan and M. Biittiker: Phys. Rev. Lett. 92 (2004) 
247901, X. Wang: Phys. Rev. A 66 (2002) 064304. 
G. Vidal, J.I. Latorre, E. Rico and A. Kitaev: Phys. Rev. Lett. 
90, (2003) 227902. 

M. Michel et al.: Europhys. J. B 34 (2003) 325. 



